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Обработка данных аэрофотосъемки 
в задачах мониторинга 
сельскохозяйственных полей 


В данной статье описаны результаты исследования применимости некоторых подходов для решения 
задачи определения заболевания сельскохозяйственных полей по данным аэрофотосъемки. Использо- 
вались следующие подходы: пороговая сегментация по оттенку и насыщенности, нечеткая сегментация, 
вычисление фрактальных сигнатур (последний подход использован к каждому из цветовых каналов 
СНИМКОВ отдельно). Также использована нечеткая сегментация по шести каналам: три цветовых канала и 
три канала, представляющие собой результат вычисления фрактальных сигнатур по цветовым каналам. 


Введение 


Дистанционные методы мониторинга сельскохозяйственных полей дают 
возможность оперативно выявить участки полей, пораженные болезнью. Выявление 
заболевания на ранних стадиях развития значительно сокращает затраты труда и 
средств. Находят развитие два основных подхода к решению задачи выявления пора- 
женных участков — спектрометрический и оптический [1], [2]. Спектральный подход 
позволяет определять многие заболевания на ранних стадиях развития. Однако, 
несмотря на это, оптические методы также существенно развиваются, так как их 
использование позволяет увеличить качество идентификации. 

Выделяются различные виды характеристик, которые можно использовать для 
идентификации заболевания [3-5]: геометрические, морфологические и цветовые, а 
также их комбинации, которые позволяют уменьшить пространство признаков, что 
упрощает идентификацию. 

Основной задачей при обработке изображений для идентификации пораженных 
болезнью участков растения или поля является сегментация, т.е. выделение однородных 
по своим цветовым или фрактальным признакам объектов изображения (группы 
пикселей) и отнесение их к тому или иному заранее определенному классу. 

Целью данной статьи является исследование применимости выбранных методов 
обработки изображений для решения задачи определения начала развития заболевания 
на поле картофеля по цветовым и фрактальным характеристикам листьев. 

В качестве исходных данных использованы снимки участка картофельного 
поля, выполненные с разных высот в период со 2 по 6 июля 2007 года. На рис. 1 
приведены примеры исходных изображений, выполненных с высоты 100, 50, 15 
метров и с близкого расстояния. 

Для определения заболевания растений на полях картофеля по данным аэрофото- 
съемки использовались следующие подходы: пороговая сегментация по оттенку и насы- 
щенности, нечеткая сегментация, вычисление фрактальных сигнатур (последний подход 
использован к каждому из цветовых каналов снимков отдельно). Также использована 
нечеткая сегментация по шести каналам: три цветовых канала и три канала, пред- 
ставляющие собой результат вычисления фрактальных сигнатур по цветовым каналам. 
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Рисунок 1 — Примеры исходных аэрофотоснимков, выполненных с разных высот 


Описание системы 


Простая обработка массива изображений может быть сопряжена с проблемой 
отсутствия связи между собственно результатами обработки и реальными картами 
полей. Для решения этой проблемы система обработки должна иметь связь с ГИС 
для более эффективного решения задачи заболевших участков сельскохозяйствен- 
ных полей и упрощения их интерпретации людьми (рис. 2). 


Модуль 
обработки 


Разбиение для 
параллельной 


Рисунок 2 — Схема системы обработки 
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Описание подходов 


Оценка визуального качества изображения 


Качество изображения определяется большим количеством технических харак- 
теристик системы: соотношением сигналлпум и статистическими характеристиками 
шума, градационными характеристиками, спектральными (цветовыми) характерис- 
тиками, интервалами дискретизации и т.д. 

Одним из параметров, которые определяют качество изображений, является 
контраст. Поскольку изображение имеет сложный сюжетный характер, это поро- 
ждает необходимость при определении его контрастности исходить из контраста 
отдельных комбинаций элементов изображения. При этом все элементы считаются 
равнозначными, и контраст каждой их пары вычисляется по формуле 


с- Л 
я > 
ВЕ: 
у 
где Г, [;- яркости элементов сюжетного изображения. 
Сюжетность изображения предполагает возможность его использования 
человеком. Поэтому при оценке контраста, как одного из параметров качества изоб- 
ражения, необходимо учитывать ряд особенностей зрительного восприятия человека. 


Рассмотрим известный эмпирический подход к оценке визуального качества 
изображения [6]. Для формирования этой оценки рассматриваются такие параметры 


изображения, как среднеарифметическое значение яркостей [, полнота исполь- 
зования градаций яркостей, резкость изображения и его обобщенный контраст. 

В целом выражение для количественной оценки визуального качества полуто- 
новых монохромных изображений записывают так: 


О=ЕКС-ТО:КО:ВО, 


где А — нормирующий коэффициент; КС — оценка контраста изображения; [О — оценка 
уровня адаптации зрительной системы; КО — полнота использования элементами 
изображения градаций яркостей; КО -— оценка резкости. 

Параметр КС дает оценку контраста изображения. Когда оценивается визуаль- 
ное качество изображения независимо от его сюжетного наполнения, тогда за КС 
принимают обобщенный контраст изображения С.еи. 


Задача сегментации 


Цель сегментации — разбиение изображения на однородные области. Однородность 
рассматривается в смысле сходства интенсивности цвета ["]. 

Математическая формулировка задачи сегментации. 

Пусть /(х,у) — дискретная функция, описывающая анализируемое изображение; 

Х- конечное множество точек плоскости, на котором определена функция /(х, у); 

5 = {55...95 — разбиение Х на К непустых связных областей 5; и 1=/,2,...,К; 

[р — предикат, заданный на множестве 3 и принимающий истинные значения 
тогда и только тогда, когда любая пара точек (х,у) из каждого подмножества 5, удов- 
летворяет некоторому критерию однородности этого подмножества. 
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Сегментацией изображения /(х,у) по предикату [р называется разбиение 
59=1{51,52,....5, удовлетворяющее условиям: 


К 
а) [15 =Х; 
1=1 
6) 5; [] 5, = ©) для любых #27; 


в) [р( ; } = «истина» для любого й; 


г) [р(5; ) 5) =«ложно» для любых #5. 


Отметим, что К областей, полученных в результате сегментации часто 
различных классов, группируются в т классов, где 2< т < (. 


Задача нечеткой классификации 


Задача нечеткой классификации ставится следующим образом: Х = (х») = 
= (хы, Хр, ..., Хр) — объекты, подлежащие классификации, где К = [1, п]; п — количество 
объектов ([8]). Количество классов равно с. Необходимо каждому элементу множества 
Х поставить в соответствие степени принадлежности к классам. Сумма степеней 
принадлежности одного элемента к классам должна равняться единице. 


Пороговая сегментация 


При пороговой сегментации выделяются участки, цветовые характеристики кото- 
рых, соответствуют пожелтевшим листьям заболевших растений, что соответствует 
следующим диапазонам значений оттенка (Вие) и насыщенности (завиганоп) [7]: 

Бое = [0,9; 1,4], заб иганоп = [45; 69]. 


Алгоритм вычисления фрактальной сигнатуры 


Суть алгоритма вычисления фрактальной сигнатуры состоит в том, что 
квантованные значения интенсивности двумерного сигнала должны располагаться 
между двумя функциями, называемыми верхней и нижней поверхностями [9]. Верх- 
няя поверхность И содержит множество точек, значения которых всегда, по крайней 
мере на один квант интенсивности, превышают интенсивность входного сигнала. 
Нижняя поверхность [, имеет значения точек, которые всегда ниже, по крайней мере 
на один квант интенсивности, входного изображения. 

Верхняя и нижняя поверхности при нулевой шкале масштабов определяется как: 


ОС, /,0) = ЦЬ 7,0) = 2,7); 
где 2(1,/) — входное изображение. 
В общем случае: 


(Е, = +1) = мах |96.1.5) - пах т, 5) 


К,те 
КЕЙ 


ГР, „5-1 =шШ |.) - пахА4,т, 2): 


п = {(К,т) [а [(,т), (6, 7151}; 


где 4- функция расстояния. 


«Штучний 1нтелект» 32008 285 


Ганченко В.В. 


Сконструированное покрытие, образованное двумя указанными функциями, 
имеет толщину 2=. Для двумерного сигнала площадь поверхности есть объем, занятый 
покрытием и деленный на величину 2=. Площадь «поверхности» интенсивности 4() в 
пределах окна наблюдения А рассчитывают вычитанием точки за точкой нижней 
«поверхности» из верхней с дальнейшим суммированием по всему окну: 


ОС, ),=) - МЕ Л,=) 
вы но 
2Е 2. 


Фрактальную размерность определяют по наклону [00° А(=) как функцию [0® 5. 


Оценка сходства изображений 


Для автоматического сравнения используют меры сходства [10]. Для сравнения 
исходного и сегментированного изображений необходимо представить последнее в 
виде псевдоцветного изображения. Для этого каждая область сегментированного 
изображения раскрашивается средним цветом соответствующей области исходного 
изображения. Сегментированные изображения сравнивают с исходным изображе- 
нием или с сегментируемым изображением идеального алгоритма, в результате чего 
для каждого алгоритма получаем числовую оценку сходства сегментируемого и 
идеального изображений. 

Функция Лью и Янг. Лью и Янг (Гла, Уап?) оценивали результат сегментации 
цветных изображений с помощью функции: 


К е 

ол УЕ 
10000%хМ) 1 /А; 
где / — сегментируемое изображение размером Мх№; К — количество областей на 
сегментируемом изображении; 47 — площадь 1-Й области; е; — средняя цветовая 
ошибка 1-й области, которая равняется сумме евклидовых расстояний между КОВ 
цветовым вектором пикселей 1-й области (на исходном изображении) и цветовым 
вектором, соответствующим области 1 на сегментированном изображении. 
Алгоритм, для которого значение функции Ё2(1) меньше, считается лучшим. 


2 


Е2(1) 


Результаты обработки 


В качестве примера на рис. 3 — 6 приведены результаты обработки для снимков, 
выполненных с высоты 50 метров. По оси абсцисс отложены дни развития 
заболевания. 


Результат оценки визуального качества изображений 


2000 | 7% —®_ Ряд1 
4 
0 | | Г. | № | № | \ 7) Ряд3 
1 2 Е 4 5 


Рисунок 3 — Оценка качества для цветовых каналов исходных изображений 
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Рисунок 5 — Оценка качества для результатов нечеткой сегментации 


Эксперимент 1- результат сегментации цветовых каналов исходных изобра- 


жений на пять сегментов; 


Эксперимент 2 — результат сегментации цветовых каналов исходных изображе- 
ний совместно с результатами вычисления соответствующих каналам фрактальных 


сигнатур на пять сегментов. 


80 - 
60 - 
40 - 
20 - 
0 


1 2 3 4 5 


—®—Ряд1 


Рисунок 6 — Оценка качества для результатов пороговой сегментации 


Результат автоматизированной оценки сходства изображений 


Результаты обработки сравниваются с исходными изображениями. 


0,06 

ВЧ Е. 

002 | 7ч —®_ Ряд2 
0 а т ы т ы -&_ РядЗ 


1 2 3 4 5 


Рисунок 7 — Оценка сходства для фрактальных сигнатур цветовых каналов 


«Штучний 1нтелект» 32008 


287 


Ганченко В.В. 


0,01 - 0,05 
0,04 - 
0,005 вы 
Ряд1 0,02 —9—Ряд1 
0,01 - 
0 71 т Т 0 
Эксперимент 1 Эксперимент 2 


Рисунок 8 — Оценка сходства для результатов нечеткой сегментации 


Эксперимент | — результат сегментации цветовых каналов исходных изоб- 
ражений на пять сегментов; 

Эксперимент 2 — результат сегментации цветовых каналов исходных изображений 
совместно с результатами вычисления соответствующих каналам фрактальных сиг- 


натур на пять сегментов. 
0,01 
—®— Ряд1 
0,005 - 


Рисунок 9 — Оценка сходства для результатов пороговой сегментации 


Выводы 


Исходя из приведенных графиков видно, что некоторые оценки позволяют 
сделать выводы относительно развития заболевания по снимкам, выполненным с раз- 
личного расстояния. Так, при оценке визуального качества исходных изображений 
отмечается уменьшение численных значений оценки по мере развития заболевания. 
Однако это наблюдается только на снимках, выполненных с расстояния 100 и 50 
метров, и проявляется наиболее ярко для красного и зеленого каналов исходных 
изображений. Схожее явление заметно и при оценке визуального качества результатов 
пороговой сегментации (для высот 100, 50 и 15 метров), а также для оценки результатов 
нечеткой сегментации для снимков, выполненных с 15 метров, причем результаты 
сходны как при сегментации цветовых каналов изображения, так и при совместной 
сегментации цветовых каналов с результатами вычисления их фрактальных сигнатур. В 
других случаях использованная оценка визуального качества изображений не дает 
видимых результатов. 
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При оценке сходства результатов обработки с исходными изображениями 
отмечается уменьшение численных значений оценки для высот 100, 50 и 15 метров. 
Оценки результатов остальных подходов к обработке изображений видимых резуль- 
татов не дают. 
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